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THE  THERMAL  ACCOMMODATION  COEFFICIENT:  A  CRITICAL  SURVEY 


Harold  Y.  Wachman* 


INTRODUCTION 

A  problem  that  currently  confronts  the  engineer  concerned  with  high- 
altitude  flight  involves  the  determination  of  drag  coefficients  for  a  vehicle  in  motion 
in  a  rarefied  gas.  The  free-molecule  drag  coefficient,  Cj^,  could,  in  principle,  be 
readily  computed  if  the  vehicle  were  to  receive  from  each  incident  gas  molecule, 
upon  collision,  all  the  momentum  that  the  molecule  has  relative  to  the  vehicle, 
i.  e. ,  if  the  incident  gas  molecules  were  brought  to  rest  (ignoring  thermal  motion) 
relative  to  the  vehicle's  surface  upon  collision  with  it.  However,  this  kind  of 
molecular  behavior  is  unusual.  It  is  more  likely  that  after  a  single  encounter 
with  a  surface,  gas  molecules  will,  on  the  average,  rebound  from  it  with  appreciable 
velocities.  Therefore,  if  the  gas  about  the  vehicle  surface  is  rarefied,  mole¬ 
cules  will  escape  the  interface  after  surrendering  to  the  surface,  on  the  average, 
only  a  fraction  of  their  incident  relative  momentum,  and  the  drag  coefficient 
will  depend  not  only  on  the  velocities  of  the  incident  molecules  but  also  on  the 
velocities  of  the  ones  rebounding. 

In  the  case  of  artificial  earth  satellites  the  choice  of  a  numerical  value  for 

the  average  speed  of  the  incident  molecules  relative  to  the  satellite  surface,  , 

presents  no  special  problem.  At  altitudes  of  less  than  1000  miles,  where  the 

atmospheric  density  is  sufficiently  high  to  produce  significant  neutral  particle 

'''Space  Sciences  Laboratory,  Missile  and  Space  Vehicle  Department,  General 
Electric  Company,  Philadelphia,  Pa. 
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drag,  the  satellite  speed  exceeds  the  mean  thermal  molecular  speed  by  a  factor 

of  6  or  more^-\  Hence  is,  nearly  enough,  equal  to  the  flight  speed  of  the 

satellite.  However,  the  mean  speed  of  the  rebounding  molecules  relative  to  the 

satellite,  V  ,  is  not  known  a  priori  and,  since  relevant  experimental  data  are 
r 

not  yet  available,  numerical  estimates  of  free  molecule  drag  coefficients  often 

lack  the  required  precision.  For  example,  data  obtained  from  observations  of 

orbital  decay  of  artificial  satellites  can  be  used  to  calculate  the  product  of 

and  the  atmospheric  density.  The  fact  that  average  values  of  the  atmospheric 

density  at  high  altitudes  derived  from  these  data  are  not  very  precise  is  due, 

(2) 

in  part,  to  the  uncertainty  in  Harris  and  Jastrow  — ^  obtained  average 

densities  using  =  2. 3  ±  0. 3,  a  value  near  the  mean  of  two  values  (2, 0  and 

2, 7)  calculated  on  the  basis  of  extreme  assumptions  regarding  the  nature  of  the 

interaction  of  gas  molecules  at  the  satellite  surface. 

Perhaps  the  simplest  way  of  relating  V  to  V.  is  by  means  of  the  thermal 

r  1 

accommodation  coefficient  a .  The  accommodation  coefficient  is  a  quantity  which 

is  used  to  express  the  efficiency  of  energy  transfer  between  a  gas  and  a  boundary; 

its  relation  to  the  ratio  V  /V. ,  and  to  C_ ,  will  be  discussed  below.  Experimental 

r  1  '  D 

studies  of  thermal  accommodation  have  been  mostly  concerned  with  gases  at,  or 
near,  room  temperature.  Therefore,  the  applicability  of  available  values  of  rt 
to  the  calculation  of  hyperthermal  free-molecule  drag  coefficients  is  quite  uncertain. 
However,  methods  which  have  been  employed  to  determine  the  thermal  accommo¬ 
dation  coefficient,  and  physical  factors  which  have  been  found  to  influence  its 
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magnitude  and  behavior  are  pertinent  to  all  future  studies.  The  literature  selected 
for  review  here  represents  efforts  which,  in  the  writer's  judgment,  have  contri¬ 
buted  most  to  an  understanding  of  Ihe  phenomenon  of  thermal  accommodation. 

THE  THERMAL  ACCOMMODATION  COEFFICIENT 

Kundt  and  Warburg'—  ,  in  an  effort  to  test  Maxwell's  hypothesis  that 
gas  viscosity  is  independent  of  density  examined  the  effect  of  changing  gas  density 
on  the  damping  of  a  vibrating  disk.  Over  a  wide  range  of  densities  the  gas 
behaved  as  predicted,  but  at  low  densities  viscosity  appeared  to  decrease  with 
decreasing  density.  Kundt  and  Warburg  suggested  that  the  interaction  bePveen 
a  low  density  gas  and  the  surface  of  a  solid  is  incomplete,  and  proposed  that  a 
gas  at  low  density  "slides"  or  "slips"  over  a  surface  which  is  transmitting  shear 
to  it.  Maxwell  showed  that  the  Kinetic  Theory  can  accoimt  for  the  slip  phenomenon 
if  in  the  treatment  of  viscous  flow  it  is  assumed  that  a  portion  f  of  every  unit 
area  of  solid  surface  "absorbs  all  the  incident  molecules,  and  afterwards  allows 
them  to  evaporate  with  velocities  corresponding  to  those  in  still  gas  at  the  tempera¬ 
ture  of  the  solid,  while  a  portion  (1-f)  perfectly  reflects  all  the  molecules  incident 

/  5  \ 

upon  it  -  Maxwell's  hypothesis  implies  that  the  energy  of  the  molecules  reflected 
by  the  portion  (1-f)  is  not  altered  during  collision,  and  suggests  a  thermal 
analogue  for  the  phenomenon  of  viscous  slip,  i.  e. ,  a  discontinuity  of  temperature 
at  an  interface  shared  by  a  low  density  gas  and  a  solid  which  are  at  different 
temperatures.  Smoluchowski'— ^  verified  the  existence  of  the  thermal  phenomenon 
when  he  measured  the  heat  conducted  by  hydrogen  from  a  heated  thermometer 
bulb  and  found  that  the  energy  interchange  between  hydrogen  and  a  glass  surface 
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is  also  incomplete.  Smoluchowski's  analysis  is  outlined  in  texts  on  the  Kinetic 
Theory  of  Gases 

In  1910  Knudsen  started  to  consider  in  detail  the  phenomenon  of  molecular 

(as  opposed  to  "molar")  heat  conduction.  He  treated  the  problem  somewhat  as 

follows  :  Consider  the  process  of  incomplete  energy  exchange  between  a  gas 

at  the  temperature  and  a  solid  at  T^;  say,  T^>  T^.  The  average  energy 

flux  through  the  surface  of  the  solid  is  the  difference  between  the  energy  brought 

up  by  an  incoming  molecular  stream  and  the  energy  carried  away  by  a  receding 

stream.  The  two  streams  may  be  considered  independently  if  the  gas  density 

is  very  low.  Accordingly,  if  the  mean  energy  per  molecule  in  the  incoming 

stream  is  e  (corresponding  to  T  )  and  N  is  the  number  of  such  molecules 
S  g 

striking  unit  surface  area  per  second,  then  the  incident  flux  of  energy  is 
Ne  =  E  .  Since  the  interaction  is  incomplete  the  mean  energy  of  a  molecule  in 

o  o 

the  receding  stream,  may  be  considered  to  be  the  same  as  the  mean  energy 

of  a  molecule  issued  from  a  body  of  gas  at  some  temperature  T  whose  value 

r 

lies  between  T  and  T  .  The  energy  carried  away  from  unit  surface  area  per 
second  is  thus  E^  =  N€^,  ^since  there  is  no  accumulation  of  gas  at  the  surface 
N  (incoming)  =  N  (receding)]  and  the  energy  flux  through  the  surface  of  the  solid 
Is  E  -  E  .  If  the  interaction  were  complete,  i.  e, ,  if  the  receding  stream  were 

O 

to  carry  the  same  mean  energy  per  molecule,  c  .  as  does  a  stream  issuing  from 

s 

a  gas  in  equilibrium  at  the  surface  temperature  T  ,  then  the  energy  flux  through 

s 

the  surface  would  be  E  -  E  where  E  =  Nc  .  These  two  energy  differences 

S  o»  S  8 

suggest  that  the  average  efficiency  of  the  energy  exchanged,  per  collision, 
between  a  gas  and  a  solid  at  their  interface  could  be  expressed  as  a  coefficient 
of  accommodation,  cc ,  defined  by, 
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E  -  E 


E  -  E 
s  g 


(1) 


Knudsen  himself  defined  a  in  terms  of  the  temperatures  of  the  molecular 
streams  as  a  measure  of  the  extent  to  which  gas  molecules  which  strike  and 
rebound  from  the  surface  of  a  solid  accommodate  themselves  to  the  temperature 


/Q  \ 

of  tile  solid.  His  formula'—'. 


a  - 


T  -  T 

T  -  T 
s  g 


(2) 


can  be  obtained  from  Equation  (1)  by  setting,  in  general,  E  =  N  (2kT)  where  k 

is  the  Boltzmann  constant  and  2kT  =  e ,  is  the  mean  translational  energy  per 

molecule  issued  from  a  body  of  gas  in  equilibrium  at  the  temperature  T. 

In  presenting  the  derivations  of  Equations  (1)  and  (2)  we  have  accepted  the 

assumption  that  the  energy  E  carried  by  the  receding  molecular  stream  is  the 

r 

same  as  the  energy  carried  by  a  stream  issued  from  a  body  of  gas  at  a  definite 

temperature  T^.  There  is  some  doubt  regarding  the  validity  of  this  assumption. 

However,  since  E^  is  a  numerical  quantity  it  can  nevertheless  be  e^qiressed  as 

E^  =  N  (2kT^)  where  is  not  a  temperature  but  a  parameter  whose  value  lies 

between  T  and  T  Alternatively,  if  a  is  determined  under  conditions  such 

s  g 

that  the  temperature  difference  =  A  T  vanishes,  the  assumption  that 

stands  for  the  temperature  of  the  receding  stream  becomes  acceptable.  Accord¬ 
ingly,  Jackson  and  Mott^~^  and  Jackson  and  Howarth^—^  employed  the  defini¬ 


tion: 


lim 

at— 0 


T  -  T 

_ g 

A  T 


(3) 
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Subsequent  studies,  which  are  cited  below,  suggest  that  within  the  range 
studied  the  magnitude  of  a  is  independent  of  ^  T.  Therefore  the  modification 
introduced  by  Equation  (3)  may  be  unnecessary. 

THE  RELATION  BETWEEN  THE  ACCOMMODATION  COEFFICIENT  AND  THE 
FREE-MOLECULE  SATELLITE  DRAG  COEFFICIENT 

Shamberg's^—^  analytical  formulation  of  the  hsqjerthermal  free-molecule 

drag  coefficient,  Cj^,  will  serve  to  illustrate  the  general  relation  between 

and  a.  The  term  "hyperthermal"  refers  to  the  fact  that  the  random,  thermal, 

motion  of  gas  molecules  is  neglected  in  the  analysis  on  the  assumption  that  the 

speed  of  the  vehicle  is  several  times  (6  or  more)  greater  than  the  mean  thermal 

speed.  The  model  consists  of  a  uniform  parallel  beam  of  molecules  that  strikes 

a  flat  elemental  surface  at  an  angle  of  incidence  9  .  and  is  then  reflected  as  a 

conical  or  wedge-shaped  beam  having  a  half-angular  width  0  and  an  axis  which 

0 

is  inclined  at  an  angle  n  ^to  the  plane  of  the  surface.  Figure  1.  Taking  all  molecular 
speeds  relative  to  the  surface  of  the  vehicle,  the  speed  of  the  molecules  in  the 
incident  parallel  beam  is  V.  (which  is  the  flight  speed  of  the  vehicle),  and  the 
speed  of  the  molecules  in  the  reflected  conical  beam  is  V^.  In  order  to  charac¬ 
terize  the  dependence  of  0  ^  upon  9^  Shamberg  chose  the  relation 

cos  =  (cos  V  2*1  (4) 


In  terms  of  these  parameters  the  hyperthermal  free-molecule-drag  coefficient 
for  all  geometries  is  expressed  by  Shamberg  as: 

r  V  1 


1  +  (5  (0  ) 

O 


f  (\;,  shape) 


(5) 
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The  second  term  inside  the  brackets  represents  the  contribution  of  the  momentum 

of  the  rebounding  molecules  to  the  drag  on  the  body.  This  term  can  be  positive 

or  negative  depending  on  whether  the  molecules  rebound  in  the  direction  of  flight 

or  in  the  opposite  direction,  respectively.  The  effective  average  direction  of 

the  rebounding  molecules  is  taken  into  account  by  the  factor  f  (v,  shape),  where 

"shape”  refers  to  the  contribution  of  the  body  shape,  and  v  comes  from  the 

reflection  law  postulated  in  Equation  (4).  The  ratio  V^/V.  is  related  by  Sham- 

berg  to  the  thermal  accommodation  coefficient  somewhat  as  follows:  if  r|  is  the 

rate  at  which  gas  molecules  collide  with  unit  surface  area,  and  m  the  mass  of 

/  2 

a  molecule  then,  in  general,  E  =  1/2  n  m  V  .  Considering  a  process  during 
which  each  molecule  makes  only  one  collision  with  the  surface,  applying  Equation 
(1)  and  using  the  notation  introduced  in  this  section  he  obtains, 


Since  V  ,  the  speed  of  the  satellite,  is  zero  in  the  chosen  reference  frame 


and 


_!r  ..  ,  1/2 

—  =  (1-a) 

i 

1/2, 


1  +  « {<^^)  ( 1  -  a)  f  {^; .  shape)  J 


(7) 


(8) 


The  following  are  some  general  inferences  drawn  in  Reference  (1)  from  a 
parametric  study  of  Equation  (8)  regarding  the  drag  coefficient  for  a  satellite  in 
the  free-molecule  regime:  (1)  The  drag  coefficient  is  generally  more  sensitive 
to  changes  in  the  nature  of  the  interaction  between  the  gas  molecules  and  the 
surface  of  the  satellite  than  it  is  to  the  geometry  of  the  satellite.  (2)  The  drag 
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coefficient  is  more  sensitive  to  changes  in  accommodation  coefficient  than  it  is 
to  changes  in  the  mean  direction  of  the  rebounding  molecules.  (3)  Without 
accurate  knowledge  regarding  the  magnitude  of  the  accommodation  coefficient 
the  detailed  composition  of  the  upper  atmosphere  may  be  neglected  in  estimates 
of  satellite  drag. 

EXPERIMENTAL  METHODS  FOR  DETERMINING  THE  ACCOMMODATION 
COEFFICIENT. 

There  are  two  important  methods  for  determining  thermal  accommodation 
coefficients.  About  two-thirds  of  the  reported  values  have  been  obtained  by 
Knudsen’s  low-pressure  method,  the  remaining  one-third  by  the  temperature- 
jump  method'  — In  the  application  of  either  method  the  experimental  cell 
usually  consists  of  a  fine  wire  filament  of  radius  r^  mounted  axially  in  a  cylin¬ 
drical  tube  of  radius  r  ,  (r  <  <r  ).  The  filament  is  maintained  at  a  constant 

^  X  A 

temperature  T^  by  passing  a  small  electric  current  through  it,  and  the  tube  with 

filament  mounted  is  immersed  in  a  bath  kept  at  a  constant  temperature  T  .  The 

w 

experimental  gas  contained  in  the  tube  is  at  a  pressure  <  0.  05  mm  when  utilizing 

the  low  pressure  method  and  roughly  in  the  range  10  -  100  mm  when  utilizing 

the  temperature-jump  method.  The  difference,  in  principle,  between  the  two 

methods  is  in  the  manner  by  which  the  energy  of  the  incident  molecular  stream 

E  ,  is  determined.  This  energy  can  be  known  in  a  simple  way  only  if  it  can  be 

g 

identified  with  the  energy  carried  by  molecules  issued  from  a  body  of  gas  at 
a  definite  temperature.  In  the  low  pressure  method  molecular  mean-free-paths 
are  large  enough  relative  to  the  relevant  dimensions  of  the  components  of  the 
cell  that  gas  molecules  transporting  energy  from  the  heated  filament  to  the  tube 
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walls  can  achieve  thermal  equilibrium  with  the  tube  walls  at  before  returning 
to  the  filament,  or  before  making  collisions  with  other  gas  molecules  or  with 
each  other.  Furthermore,  as  a  result  of  the  cyllndrically  S3mimetric  cell 
geometry,  and  the  fact  that  r  ■<s<  r  ,  the  number  of  molecules  which  at 

X  A 

any  given  time  may  find  themselves  in  the  receding  stream  moving  from  the  fila¬ 
ment  to  the  tube  walls  with  excess  energy  acquired  at  the  filament  is  a  very 
small  fraction  of  the  total  number  of  molecules  in  the  cell,  (This  fraction  can 
be  shown  to  be,  roughly  r^ 

Consequently  the  entire  body  of  gas  in  the  tube  may  be  considered  to  be  at  the 
temperature,  T^,  of  the  tube  walls,  thus  =  T^.  At  the  gas  pressures  utilized 
for  the  temperature -jump  method  the  conduction  of  heat  from  filament  to  tube 
walls  induces  temperature  gradients  in  the  gas  and  the  assumption  is 

not  valid.  Instead  the  method  relies  on  a  relation  between  the  accommodation  coef¬ 
ficient  and  the  temperature  jump  distance;  the  latter  quantity  is  determined  from 
experimental  measurements. 

The  detailed  experimental  procedure  employed  varies  from  laboratory  to 
laboratory,  however,  the  following  is  common  in  many  investigations  with  either 
method.  After  the  experimental  tube  is  joined  to  a  vacuum  system  it  is  sur¬ 
rounded  by  a  bath  kept  at  a  constant  temperature,  T^,  and  evacuated.  The  gas 
under  investigation  is  introduced  into  the  tube  at  the  desired  pressure  and  the 
filament  is  heated  electrically.  The  current  required  to  maintain  the  filament 

at  an  average  temperature,  T  ,  and  the  potential  drop  across  the  filament  are 

s 

measured,  and  these  data  are  used  to  calculate  the  power  input  to  the  filament 
Wj,  and  its  resistance.  The  average  temperature  of  the  filament  is  calculated 


when  the  mean  free  path  is^  ^r^  ~ 
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accurately  from  its  resistance  and  the  temperature  coefficient  of  resistance  of 

the  material.  W  ,  the  total  power  required  to  maintain  the  filament  at  T  ,  is 
f  s 

equal  to  the  power  which  the  filament  loses  by  radiation,  end  conduction  to  the 

supports,  and  gas  conduction.  The  power  lost  by  radiation  and  end  conduction 

is,  approximately,  W  ,  the  power  which  the  filament  kept  at  T  would  lose  in 

V  s 

vacuum.  Therefore,  W  is  also  measured  and  its  value  subtracted  from  W,  to 

V  f 

obtain  ,  the  net  steady  state  power  loss  from  unit  filament  area  to  gas, 

A.  The  Low-Pressure  Method 

The  accommodation  coefficient  is  computed  using  an  expression  introduced 

(14) 

as  a  definition  for  a  by  Blodgett  and  Langmuir —  ; 


0.  = 


W 

_£ 

W, 


KT 


(9) 


where  =  E  -  E  is  the  free  molecule  heat  conductivity  calculated  from  the 
KT  s  g 

Kinetic  Theory  as  follows:  As  shown  above,  E  in  general  may  be  expressed  as, 
E  =  2  N  k  T  (10) 

The  Kinetic  Theory  expression  for,  N,  the  number  of  molecules  striking  (or 
rebounding  from)  unit  filament  area  per  second  is 

N  =  nc/4  (11) 

where  n  is  the  molecular  number  density  and  c  the  average  molecular  velocity, 
(both  quantities  are  referred  to  a  gas  at  T^).  Converting  to  molar  quantities 


with  R  as  the  molar  gas  constant,  A  Avogadro's  number,  M  the  molecular 

weight  of  the  gas,  and  using  the  relations  n  =  AP/RT  and  c  =  (8RT  /nM) 

g  K 


1/2 


becomes 


/ 


2  TT  M  R  T 


g 


(12) 
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The  quantity  2R  is  the  translational  molar  heat  capacity  of  a  molecular  stream 
(it  is  greater  than  the  corresponding  quantity  for  still  gas  by  R/  2) .  A  more 
general  expression  is  obtained  by  replacing  2R  by  C  +  R  /  2  where  C  is  the 

V  V 

molar  heat  capacity  at  constant  volume.  Incorporating  this  into  Equation  (12) 

( 15 ) 

and  using  A  T  =  -T^  and  rearranging  results  in  ' — 


Values  of  A  ,  the  free-molecule  heat  conductivity  per  unit  AT  unit  pressure 

and  unit  filament  area,  for  several  gases  at  30°C  are  listed  in  Ref,  (15). 

B.  The  Temperature-Jump  Method 

(  7  ) 

Kennard  ~  has  outlined  the  temperature-jump  method  in  the  form  applied 

by  a  majority  of  modern  investigators.  The  method  utilizes  measurements  of 

the  thermal  conductivity  of  a  gas  at  a  series  of  pressures  which  lie  in  a  range 

within  which  the  conductivity  should  be  independent  of  pressure.  However,  the 

measured  conductivity  appears  to  diminish  with  decreasing  pressure  due  to  the 

pressure  dependence  of  the  temperature  discontinuity  at  the  gas  boundaries.  In 

computing  accommodation  coefficients  from  temperature-jump  data  obtained  in 

cells  of  cylindrical  symmetry  the  temperature  discontinuity  at  the  tube  walls 

may  be  neglected' —  .  The  accommodation  coefficient  is  related  to  the  tempera- 

( 7 )  ( 13  ) 

ture-jump  distance  at  the  filament,  g  ,  by 

s 


2-a 


(14) 


where  K_  is  the  (true)  thermal  conductivity  of  the  gas  at  the  temperature  T  , 
T  e 

e 

and  T  (Kennard ’s  T  )  is  the  temperature  that  would  obtain  in  the  gas  at  the 

6  K 
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filament  surface  if  the  temperature  gradient  (dT/dr)^  which  exists  in  the  body 
of  gas  at  a  distance  r  from  tube  axis  were  to  continue  in  regular  functional 
relation  to  r  right  up  to  the  filament  itself  (r  =  r^).  The  temperature-jump 
distance  g  is  obtained  as  follows:  The  temperature  discontinuity  T  -  T  is 

S  S  6 

related  to  g  by  the  Poisson  temperature-jump  equation: 
s 


T  -  T  =  -  g 
s  e  s 


f— ) 

\d  r  /r  = 


(15) 


y  20  \ 

Smoluchowski'  —  ^  showed  that  g  is  inversely  proportional  to  the  gas  pressure; 

s 

(17)  (13) 

thus  the  product  g  P  in  Equation  (14)  is  often  replaced  ' — '  ' — '  by  g'  which  is 
s  s 

independent  of  pressure.  The  heat  conducted  in  a  cylindrical  cell  per  unit  fila¬ 
ment  area  per  second,  W  ,  is  obtained  by  combining  the  Poisson  e^qiression 

S 

with  the  Fourier  heat  conduction  equation  to  give'-'  ' — 


W 

_-g. 

AT 


K 


^  In  (r/r^)  +  g^/P 


(16) 


where  K  is  the  thermal  conductivity  at  some  temperature  between  T  and  T  . 

8  W 

Experimental  measurements  determine  T  and  T  ,  and  W  as  a  function  of  P, 

s  w’  g 

In  the  modern  work  which  begins  with  Dickins'  — '  l/v  is  plotted  against  1/P  to 

obtain  a  straight  line  according  to  Equation  (14)  of  the  form  y  =  mx  +  b,  with 

slope  m  =  g'  /K  and  intercept  b  =  (r,  /K )  ln(r-/r,).  The  value  of  K,  determined 
s  1  Z  1 

from  the  intercept,  when  mutliplied  by  the  value  of  the  slope  gives  g*  . 

s 

A  computation  of  a  by  Equation  (14)  also  requires  an  estimate  of  T  and 

e 

( 13 )  — 

K  .  Thomas  and  Golike'  — '  proceeded  as  follows:  Setting  g  P  =  m  K  they 
1.  s 

e 

obtain: 
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2 -a  K 

a  ^  K 


Te 


V 

R 


=  g  m 


(17) 


Within  a  narrow  temperature  range  K,  in  general,  is  proportional  to  T  where  a 


is  a  constant  whose  value  depends  on  the  gas.  The  temperature  T_  at  which 
_  K 

K  is  the  conductivity  is  found  to  be  (for  He,  Ne  and  CO  )  very  close  to  the  mid- 

temperature  of  T  and  T  ;  i.  e.  T— =T  +A  T  /2  where  AT  =T  -T  . 

ew  Kwe  eew 

With  these  relations  8  in  Equation  (15)  becomes: 

/  AT^, 

C 


8  = 


IT  + 
\  w 


1/2 


T  +  AT 
w  e 


a  +  1/2 


V  ^  1 


R 


2  R 
M 


(18) 


( 13 ) 

The  quantity  AT  is  approximated  as  AT' — This  approximation  is 
e 

better  for  Ne  than  for  He  and  best  for  CO^. 

C.  Comparison  of  the  Methods 

The  low-pressure  method  is  based  on  a  theory  which  becomes  quite  accurate 
at  very  low  pressures.  Experimental  results  which  are  cited  in  subsequent  sec¬ 
tions  demonstrate  that  when  the  experimental  details  are  properly  attended  the 
low  pressure  method  gives  reproducible  results  which  are  consistent  from 
laboratory  to  laboratory. 

The  temperature  jump  method  is  not  as  well  founded.  Theory  relating  a 

to  g  is  subject  to  some  uncertainty  and  Situation  (14)  is  only  an  approximate 
s 

(7)  (17) 

relation  between  these  quantities  -  .  Grilly,  Taylor  and  Johnston —  do  not 

consider  a  obtained  by  Equation  (14)  to  be  the  same  as  the  Knudsen  low-pressure 

accommodation  coefficient.  They  write:  "A  theoretical  relation  between  g*  and 

s 

Knudsen' s  accommodation  coefficient  has  not  as  yet  been  obtained  and  probably 
no  simple  relation  between  them  exists,  because  of  the  complex  nature  of  the 
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actual  stream  of  incident  molecules^—)  They  define  the  quantity  a  of  Equation 
(14)  not  as  the  accommodation  coefficient  but  a  slightly  different  coefficient  which 
they  name  the  "Maxwell-Knudsen”  accommodation  coefficient.  There  is,  in  fact, 
very  poor  agreement  between  the  two  methods  on  the  basis  of  comparison  of 
published  values  of  accommodation  coefficients'  —  'and  until  recently  the  ability  of 
the  two  methods  to  give  the  same  a  value  was  questionable.  In  1954  Thomas 
and  Golike'  — '  published  the  results  of  a  comparative  study  of  the  accommodation 

coefficients  of  He,  Ne  and  CO  on  platinum  by  the  two  methods.  Most  of  the 

^  ^  •  1  .  „,(18)  (19)  (20)  (21)  (17)  (13  )  _ 

modern  temperature-jump  work  is  on  Pt —  —  —  —  —  —  .  The 

comparison  was  made  using  the  same  tube  and  vacuum  system,  with  gases  from 

the  same  reservoir  and  purification  system.  Within  the  limitations  imposed  by 

the  requirements  of  the  methods,  experimental  procedures  and  conditions  were 

maintained  the  same.  Agreement  with  He  and  CO  was  very  close.  For  Ne 

A 

the  results  were  ambiguous,  but  the  discrepancies  (which  were  of  the  order  of 
25%)  were  associated  with  consistent  failure  to  reproduce  the  same  surface  con¬ 
ditions  rather  than  with  failure  of  the  temperature-jump  method. 

EXPERIMENTAL  PROCEDURES 
A,  Determination  of  the  Power  Loss  in  Vacuum 

Because  the  temperature  distribution  along  a  filament  is  altered  by  gas 
f  15) 

conduction' —  the  measured  power  loss  from  a  filament  in  vacuum,  W  ,  is 

V 

not  exactly  equal  to  the  radiation  and  end-condution  losses  from  a  filament 
surrounded  by  gas.  When  A  T  is  small  the  altered  distribution  induces  only  a 
small  change  in  the  radiation  loss,  the  principle  change  being  in  end  losses 
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arising  from  new  temperature  gradients  at  the  filament  ends.  For  example, 
for  a  5  mil  potassium  covered  tungsten  filament  20  cm  in  length  operated  at  A  T  — 
10°C,  c:r300°K,  neglecting  to  correct  for  the  change  in  end  losses  results 

in  an  error  of  +  14%  in  the  accommodation  coefficient  of  helium' — ' .  Attempts 
to  minimize  or  correct  the  error  have,  in  general,  been  successful.  One 
approach  is  to  use  power  loss  data  obtained  from  the  central  segment  of  the  fila¬ 
ment  only  where  changes  in  gradients  with  gas  conductions  are  diminished.  For 
example,  Thomas  and  Olmer' — '  and  also  Thomas  and  Brown' — '  attached  potential 

leads  at  a  distance  from  each  end  of  the  filament  and  measured  the  power  input 

(24) 

to  the  central  segment  alone.  Thomas  and  Schofield  — '  using  a  single  potential 

lead  subtracted  the  appropriate  measurements  taken  on  the  shorter  section  from 

those  obtained  on  the  longer  section  to  obtain  net  values  which  apply  to  the  middle 

( 13 ) 

section.  A  more  elegant  approach  was  taken  by  Thomas  and  his  co-workers' —  , 

—  ,  — '  when,  with  auxiliary  heaters,  they  supplied  heat  at  a  controlled  rate 

to  the  filament  ends  and  thus  obtained  a  flat  temperature  distribution  over  the 
entire  length  of  the  filament.  In  addition  to  these,  R.  E.  Brown  developed  an 
analytical  method  for  filament  end-loss  correction  which  is  reproduced  in  Appen- 
dix  A  to  a  dissertation  by  H.  L.  Petersen'  — ' , 

B.  Preparation  of  the  Experimental  System 

The  magnitude  of  the  accommodation  coefficient  may  be  markedly  affected 
by  gaseous  layers  adsorbed  on  the  filament  surface  formed  sometimes  quite 
slowly  from  traces  of  impurities  in  the  vacuum  system  or  from  the  gas  under 
consideration.  In  order  to  subject  these  effects  to  reproducible  control  rigorous 
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(25) 

experimental  procedures  must  be  adopted.  J.  K.  Roberts  ' — '  was  the  first  to 
recognize  this  need  and  his  pioneering  work  which  began  in  1930  is  the  basis  of 
all  modern  investigations.  However,  many  workers  who  followed  Roberts  have 
not  adopted  the  required  procedures  and  made  no  serious  attempt  to  create  or 
maintain  solid  surfaces  of  known  composition.  As  a  result,  many  published 
values  of  accommodation  coefficients  cannot  be  associated  with  the  interaction 
of  known  gases  on  known  surfaces,  and  values  for  presumably  the  same  gas  and 
surface  may  be  widely  different.  Compare,  for  example,  the  accommodation 
coefficient  values  of  helium  on  tungsten  of  0. 0165  reported  by  Thomas  and  Scho- 
field' — ‘  and  about  0.  5  obtained'  —  ^  on  uncleaned  tungsten.  The  lower  value 
was  obtained  after  the  vacuum  system,  experimental  tube,  and  filament  were 
all  subjected  to  careful  outgassing  and  the  tube  walls  covered  with  an  evaporated 
aluminum  film  to  act  as  getter.  This  value  is  very  probably  characteristic  of 
the  accommodation  coefficient  of  helium  on  tungsten.  The  value  0.  5,  however, 
cannot  be  associated  with  a  known  surface. 

It  may  be  that  a  clean  metal  surface  can  be  maintained  only  in  the  presence 
of  a  few  gases,  helium  and  neon  for  example.  However,  even  when  the  gas  under 
investigation  itself  adsorbs  on  the  surface,  the  surface  should  be  initially  clean 

since  the  magnitude  of  n:  may  depend  on  the  type  and  quantity  of  adsorbed  gas. 

( 27 ) 

Wachman  —  found,  for  example,  that  the  value  of  a  for  hydrogen  on  tungsten 
initially  covered  with  oxygen  is  appreciably  different  than  a  for  hydrogen  meas¬ 
ured  on  initially  clean  tungsten;  hydrogen  was  presumably  adsorbed  on  both  bare 
and  oxygen  covered  sui’faces  during  the  measurements. 

The  problem  of  producing  and  maintaining  a  clean  surface  is  not  an  easy 
one  since  such  a  surface  is  very  reactive  and  rapidly  acquires  an  adsorbed  layer 
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even  when  exposed  to  a  good  vacuum  (say  10  mm  of  Hg).  Therefore  the 

experimental  system  must  be  treated  to  reduce  the  quantity  of  residual  gaseous 

impurities  to  a  level  at  which  the  filament  surface,  once  cleaned,  remains  essen- 

(24) 

tially  free  of  impurities  during  the  course  of  a  run.  Thomas  and  Schofield 
preceded  as  follows:  "The  system  was  all  Pyrex  with  two  mercury  condensation 
pumps  in  tandem,  with  no  stopcocks  in  the  highvacuumline.but  with  several  special 
grind  mercury  seal  stopcocks  in  the  gas  supply  and  purifying  system,  which  was 
separated  from  the  main  system  at  all  times  by  a  mercury  cut-off  except  during 
preliminary  evacuation  and  admission  of  gas. "  They  describe  the  outgassing 
procedure  as  follows  "Ovens  are  then  placed  about  the  tube,  charcoal  trap 

on  the  gas  supply  manifold  and  liquid  air  trap.  These  are  baked  at  about  400°C 
for  several  hours,  open  to  the  pumps.  With  the  ovens  hot,  the  rest  of  the  system 
is  then  heated  repeatedly  with  a  large  brush  gas -air  flame  for  two  to  three  hours, 
almost  to  the  softening  point  of  the  glass,  with  special  attention  devoted  to  the 
region  between  the  liquid  air  trap  and  the  tube.  After  cessation  of  flaming,  the 
ovens  are  turned  off,  liquid  air  is  put  on  the  single  U-tube  on  the  mercury  side 
of  the  trap,  and  the  filament  is  heated  to  flashing  temperature." 

Apart  from  the  problem  of  preparing  the  experimental  system  the  filament 
surface  itself  must  be  cleaned.  Roberts  outgassed  tungsten  filaments  by 
heating  to  2000°C.  However,  this  procedure  can  be  applied  only  to  a  few  metals. 
The  oxide  coating  which  covers  most  metal  surfaces  is,  in  most  cases,  not 
volatile  below  the  metal’s  melting  point.  Eggleton  and  Tompkins'  —  ^  success¬ 
fully  reduced  the  oxide  coating  on  an  iron  filament  by  heating  in  hydrogen,  but 
found  that  the  process  left  a  hydrogen  film  on  the  surface  which  could  not  be 
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desorbed  by  prolonged  heating  in  vacuum.  To  obtain  a  clean  iron  surface  they 

bombarded  the  filament  with  ions  of  nitrogen  and  argon  after  the  hydrogen  reduc- 

tion.  Faust'  — '  working  with  aluminum  attempted  to  form  a  clean  surface  by 

evaporating  aluminum  onto  platinum  and  tungsten  filaments,  Brown^~^  applied 

(22) 

this  technique  to  form  clean  beryllium  surfaces,  and  Petersen'  — '  to  the  determi¬ 
nation  of  accommodation  coefficients  on  clean  sodium  and  potassium. 

RESULTS  OF  EXPERIMENT 

The  work  described  in  subsequent  sections  was  by  the  low-pressure  method 
unless  otherwise  specified. 

A.  Effects  of  the  Surface  on  the  Thermal  Accommodation  Coefficient 

Experimental  investigations  in  which  serious  attempts  were  made  to  work 

(25) 

with  surfaces  of  known  composition  begins  with  Roberts  ' — ' .  He  measured 
accommodation  coefficients  of  helium^~^  and  neon^~^  on  cleaned  tungsten  and 
obtained  values  which  were  roughly  one-fifth  the  lowest  previously  reported. 

Before  each  run  Roberts  outgassed  the  surface  by  heating  and  during  each  run  the 
helium  and  neon  were  continuously  circulated  through  charcoal  at  the  temperature 
of  liquid  air  to  remove  adsorbable  impurities  coming  from  the  glass  walls  of  the 
apparatus.  Mann' —  ,  working  with  helium,  and  Mann  and  Newell'  — '  with 
hydrogen  and  deuterium  on  platinum  also  attempted  to  clean  the  surface.  How¬ 
ever,  they  did  not  purify  the  gases  they  used  and  contaminants  carried  by  these 
adsorbed  rapidly  on  the  filament. 

Roberts'  values  for  the  accommodation  coefficient  of  helium  on  tungsten, 
n  (He  -  W)j  ,  at  300  K  lie  between  0.  05  and  0.  06,  while  values  on  untreated 
tungsten  surfaces  lie  about  0.  3.  The  fact  that  the  accommodation  coefficient 
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/  32 ) 

increases  with  adsorption  on  a  surface  was  confirmed  by  Mann'  ' ,  but  questioned 
(34) 

by  Bremner'  — '  who  ascribed  the  increase  of  a  with  time,  observed  by  Roberts, 

to  persistence  of  thermal  effects  of  flashing.  However,  Roberts'  observations 

( 35 ) 

of  the  behavior  of  a  were  reconfirmed  by  Eggleton,  Tompkins  and  Wanford'  — 

and  by  J.  W.  Faust'  —  \  Faust,  working  with  inert  gases  on  evaporated  aluminum 

surface  found  that  low  a  values  can  be  obtained  on  platinum  and  tungsten  filaments 

kept  hot  during  the  evaporation  of  aluminum  from  other  filaments  to  the  tube 

walls.  Apparently  the  aluminum  coating  on  the  tube  walls  was  acting  as  a  getter, 

thus  maintaining  a  lower  pressure  of  impurities  in  the  tube  than  could  be  obtained 

(24) 

with  the  vacuum  pumps  alone.  Thomas  and  Schofield  improved  this  technique 
and  applied  it  to  measurements  of  a  (He  -  W)  from  83^  to  303 °K.  In  all  previous 
attemtps  to  obtain  clean-surface  a  values,  the  measured  a.  increased  rapidly 
with  time  after  flashing  the  filament.  The  increase  with  time  was  caused  by 
adsorption  of  gaseous  impurities  not  entirely  eliminated  from  the  experimental 
system.  To  obtain  clean  surface  values,  Roberts  plotted  rt  vs.  time  after  flashing 
and  extrapolated  to  zero  time.  In  the  work  of  Thomas  and  Schofield  the  change 
with  time  was  almost  entirely  eliminated,  and  the  measured  values  of  a  (He  -  W), 
about  0.0165,  were  lower  than  any  reported  before.  Silvernail'  — '  virtually 
eliminated  the  change  with  time  using  a  mischmetal  getter  in  an  investigation  of 
a  (He,  Ne,  A  -  W)  from  77^  to  303‘^K. 

The  use  of  a  getter  on  the  tube  walls  appears  an  excellent  way  to  supple¬ 
ment  the  general  precautions  required  to  obtain  n.  values  of  inert  gases  on  clean 
surfaces.  However,  the  getter  technique  for  removing  adsorbable  impurities 
is  impractical  when  the  gas  whose  is  being  measured  itself  adsorbs  on  the 
getter.  When  this  is  the  case,  satisfactory  vacuum  conditions  may  be  obtained  by 
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prolonged  and  careful  outgassing  [see,  for  example,  Footnote  (7)  Reference 
and  also  Reference  (^)]  . 

In  1935  Roberts'  — '  suggested  that  the  accommodation  coefficient  of  an 

inert  gas  such  as  neon  may  be  used  as  an  indicator  for  adsorption  of  other  gases. 

He  proceeded  to  investigate  the  adsorption  of  hydrogen' —  ,  oxygen'  — '  and 
(39) 

nitrogen  on  initially  clean  tungsten  by  observing  change  in  the  accommodation 

coefficient  of  neon.  Figure  (2)  is  an  example  of  the  data  he  obtained.  The  Figure 

shows  plots  of  a  (Ne)  at  two  temperatures  vs  time  (after  flashing)  and  illustrates 

the  change  in  a  (Ne)  caused  by  the  admission  of  small  quantities  of  hydrogen  to 

the  experimental  tube.  A  summary  of  Roberts'  work  is  contained  in  his  book^  — \ 

(41) 

and  in  the  revised  edition  by  Miller  — Others  following  Roberts  have  applied 
the  accommodation  coefficient  method  with  neon  to  study  adsorption  on  initially 

bare  tungsten.  Thus  Van  Cleave'  — '  studied  adsorption  of  nitrogen,  and  Morrison 

( 43  )  ( 44 )  (  27  ) 

and  Roberts  —  —  investigated  adsorption  of  oxygen.  Wachman'  —  'using  helium 

as  an  indicator  gas  investigated  adsorpiton  of  H  ,  D„,  O  ,  N  ,  CO  ,  CH  ,  C  H  and 

/  28  ) 

C^Hg.  Eggleton  and  Tompkins'  — '  applied  the  method  to  an  investigation  of 
adsoi'ption  of  H  ,  O  and  N  or  initially  clean  iron,  and  Petersen —  to  a  study 
of  adsorption  of  potassium  on  tungsten.  The  remarkable  dependence  of  the  accom¬ 
modation  coefficient  on  the  species  of  gas  adsorbed  on  a  surface  is  illustrated  by  the 

difference  between  a  (He)  values  found  on  hydrogen  and  deuterium  covered  tungsten. 
(27) 

Wachman  — 'starting  with  a  value  of  0.0200  for  a(He)  on  bare  tungsten  measured 
0.0407  for  n(He)  on  hydrogen  covered  tungsten  as  against  0.0458  on  deuterium 
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covered  tungsten.  In  addition  a(H  )  and  a(D  )  were  also  found  to  have 
different  values  on  surfaces  of  hydrogen  and  deuterium,  respectively,  adsorbed 
on  initially  bare  tungsten. 

The  magnitude  of  the  accommodation  coefficient  depends  not  only  on 
the  condition  but  also  on  the  type  of  surface.  Tungsten  is  preferred  in  many 
investigations  because  it  can  be  obtained  in  fine  filament  form  and  is  readily 
outgassed  by  heating.  These  properties  of  tungsten  have,  in  part,  also  been 
responsible  for  its  extensive  use  in  studies  of  various  surface  phenomena 
by  other  methods  since  the  earliest  work  of  Langmuir,  ' — '  ' —  .  This 
coincidence  is  convenient  in  interpreting  results  of  experimental  and  analytical 
investigations  on  accommodation  coefficients.  There  is  at  the  same  time, 
and  partially  for  the  same  reasons,  an  extensive  bibliography  on  platinum. 
Most  of  the  temperature  jump  work,  for  example,  is  with  platinum  (see 
above).  However,  in  most  studies  with  platinum  little  effort  was  made  to 

/  23  \ 

rid  the  system  of  adsorbable  impurities.  The  work  of  Thomas  and  Brown' —  , 
though  not  achieving  the  excellent  conditions  obtained  in  later  work  with  tungsten 
' — ' —  ,  may  be  cited  as  an  exception.  The  measurements  were  with 
He,  Ne,  A,  Kr,  H  ,  N  ,  O  ,  CO,  and  Hg  on  platinum, 

Z  A  z 

Table  I  presents  a  list  of  accommodation  coefficient  values  of  helium 
and  neon  on  different  surfaces  at  about  room  temperature.  The  table  Includes 
mention  of  the  method  used  to  prepare  each  final  surface  and  the  reference. 
The  values  in  the  table  are,  in  our  opinion,  among  the  most  reliable  ones 
available. 
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Table  I 


Accommodation  Coefficients  of  Helium  and  Neon  on 
Clean  Metal  Surfaces 

Preparation  of  Accommodation 


Gas 

Metal 

Surface 

Coefficient 

Ref. 

Helium 

Tungsten 

flashing 

0.017 

(24),  (36) 

Neon 

Tungsten 

flashing 

0.042 

(36) 

Neon 

Iron 

Reduction  with  H^ 

0.056 

(28) 

and  ion  bombardment 

Helium 

Platinum 

flashing 

0.038 

cited  in  (10) 

Helium 

Aluminum 

Evaporation  on  platinum 

0.073 

(29) 

Neon 

Aluminum 

Evaporation  on  platinum 

0.159 

(29) 

Helium 

Beryllium 

Evaporation  on  tungsten 

0. 145 

(10) 

Neon 

Beryllium 

Evaporation  on  tungsten 

0.315 

(10) 

Helium 

Potassium 

Evaporation  on  tungsten 

.083 

(22) 

Neon 

Potassium 

Evaporation  on  tungsten 

0.199 

(22) 

Helium 

Sodium 

Evaporation  on  tungsten 

.090 

(22) 

Neon 

Sodium 

Evaporation  on  tungsten 

.198 

(22) 

The  following  is  a  brief  list  of  references  (n,  v. )  to  measurements  of  a  on 
surfaces  which  may  cr  may  not  have  been  characteristic  of  the  bare  material: 
a  (He  -  Ni),  ^  — a  (Hydrocarbon  molecules  -  Ni),  ^  n  (H^  -  Fe)/— 

^(He,  Ne,  H  ,  N  -  glass)  ^  ^ ;  a  (air  -  several  engineering  surfaces),  ^  ^ 

B.  The  Temperature  Dependence  of  the  Accommodation  Coefficient 

The  simplest  way  of  studying  the  effect  of  temperature  on  the  accommodation 
coefficient  is  with  a  clean  surface.  A  gas  covered  surface  usually  changes  in 
composition  with  changing  surface  temperature  and  may  also  change  with  changing 
gas  temperature.  Therefore,  observed  changes  in  a,  with  temperature,  on  a  gas 
covered  surface  may  be  due  to  changes  in  surface  composition.  An  interesting 
collection  of  data  is  presented  by  Thomas  and  Olmer'  — '  who  assembled,  in  a 
single  graph,  measured  values  of  t  vs.  T  fi*om  several  sources  [see  also  page 
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129  Reference  (8  )]  .  All  the  data  cited,  with  one  exception,  show  a  decreasing 
with  increasing  temperature  below  500*^K.  The  exception  is  a  plot  of  Roberts'  ^ 
values  of  a  (He)  on  cleaned  tungsten  which  shows  a  increasing  with  temperature. 

It  appears  a  good  assumption  that  in  all  the  work  cited  [except  Reference  (^)]  the 
surface  at  each  temperature  was  saturated  with  adsorbed  gas.  Since  the  (saturation) 
quantity  of  gas  on  a  surface  usually  decreases  with  increasing  temperature  the 
corresponding  a  values  would  also  decrease.  The  surface  in  Reference  (30  ) . 
however,  was  partially  clean,  with  surface  population  certainly  below  the  satura¬ 
tion  value.  In  this  case,  an  increase  in  T  (say,  at  constant  A  T)  may  have  increased 

S 

the  rate  at  which  contaminant  gas  from  the  tube  walls  (T  =  T  )  arrived  at  the 

g  w 

partially  covered  filament  surface  and  a  increased  because  the  surface  con¬ 
tamination  increased.  Evidently  none  of  these  data  offers  much  insight  into 
the  temperature  dependence  of  the  accommodation  coefficient. 

The  accommodation  coefficient  may  depend  on  both  T^  and  T^  (or,  more 

conveniently  on  T  and  A  T).  Among  the  most  reliable  investigations  of  this 
g 

dependence  are  the  studies  described  in  References  (24)  and  (^)  of  accomoda¬ 
tion  coefficients  of  inert  gases  on  clean  tungsten.  Figure  3  is  a  plot  of  a  (He)  vs. 


bath  temperature,  T  ,  (T  =  T  )  according  to  (24),  (curve  B),  and  according 
w  g  w 

to  (30),  (curve  A).  Curve  B  values  were  with  A  T  nearly  constant  at  about  18°; 
the  change  of  a  (He)  with  T^  is  sensible  but  small  with  values  0.0151  at-190°C, 

0.  0162  at  -135°  and  -80°C,  0. 0172  at  -30°C  and  0.  0167  at  30°C.  Silvernail^  — ^ 
(not  shown)  found  a(He)  0.016  at  30°C  and  0.  012  at  -196°C.  The  change  in  a(Ne) 


with  T^^  is  somewhat  more  marked  in  the  same  temperature  range 
exhibiting  a  slight  minimum  at  -180°C.  Figure  4,  ,  is  a  plot  of  a  (He)  values 


23 


against  time  after  the  filament  was  outgassed  by  flashing.  The  results  of  two 

runs  are  plotted  with  T  =  -190°C  and  T  changing  at  intervals  shown  on  the 

w  s 

plots;  thus  the  plots  show  the  effect  of  changing  A  T  (at  constant  T^)  on  the  mag¬ 
nitude  of  a  (He).  These  results  combined  with  other  from  Reference  (^)  indi¬ 
cate  that  for  He  and  Ne  on  clean  tungsten  a  is  nearly  constant  with  changing  A  T, 

and  changes  with  T  at  constant  A  T. 
w 

C.  The  Effect  of  Gas  Pressure  on  the  Accommodation  Coefficient 

Since  surface  composition  may  be  altered  by  changing  ambient  gas  pres¬ 
sure  the  accommodation  coefficient  may  appear  to  change  with  pressure.  How¬ 
ever,  the  accommodation  coefficient  defined  as  a  measure  of  the  efficiency  of 
the  energy  exchanged,  per  collision,  is  by  definition  independent  of  pressure. 

This  behavior  with  pressure  is  indeed  found  when  adequate  care  is  taken  to  main¬ 
tain  constant  surface  composition  in  the  course  of  experiment. 

Amdur,  Jones  and  Pearlman' — '  reported  that  the  accommodation  coeffi¬ 
cients  of  several  gases  including  He,  Ne,  A,  Kr  and  Xe  on  platinum  at  room 

/  54 ) 

temperature  changed  with  pressure  below  0. 1  mm  Hg.  Amdur  —  attributed 

this  pressure  dependence  to  simple  Langmuir  adsorption  on  the  surface,  of  the 

gases  whose  a  was  being  measured.  Morrison'  — '  accepted  this  interpretation 

for  H  ,  D  ,  N  ,  CO  and  O  but  offered  an  alternative  one  for  the  pressure  depen- 

dence  observed  with  the  inert  gases.  These  were  about  99. 9%  pure  and  Morrison 

suggested  that  the  impurities  carried  by  these  were  of  sufficient  quantities  to 

( 23 ) 

cause  the  observed  affect.  Thomas  and  Brown'  — '  measured  the  accommodation 

( 53 ) 

coefficients  of  the  gases  used  by  Amdur  et.al. ' — 'on  platinum  as  a  function  of 
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pressure  but  attempted  to  maintain  the  composition  of  the  surface  constant  by 
designing  the  experiments  to  minimize  and  control  the  quantities  of  adsorbable 
impurities  in  the  experimental  system.  Some  of  their  results  are  shown  as  plots 
in  Figures  5  and  6.  These  show  that  the  accommodation  coefficients  of  the  nine 
gases  are  essentially  constant  with  pressure  in  the  pressure  range  where  the 
low-pressure  method  is  applicable.  The  lower  values  at  higher  pressure  in 
Figure  5  are  not  real  values  of  a,  but  were  included  to  show  for  each  gas  (and 
experimental  tube)  the  range  of  application  of  the  low-pressure  method. 

Harnett criticized  the  manner  by  which  the  low-pressure  method  has 

been  applied  by  a  majority  of  the  investigators.  Rewrites:  "To  summarize,  an 

accurate  measurement  of  the  accommodation  coefficient  using  a  thin  heated  wire 

in  a  cylindrical  enclosure  must  satisfy  the  following  requirement:  ....  (2)  The 

pressure  level  within  the  enclosure  should  be  low  enough  to  ensure  that  the  mean- 

free-path  (  X)  is  at  least  three  times  the  radius  (r^)  [our  r^]  of  the  outer  cylinder, 

X  >3  r^.  (3)  The  pressure  level  within  the  enclosure  should  be  high  enough  to 

ensure  that  the  radiation  heat  transfer  (Q  )  does  not  overwhelm  the  conduction 

R 

contribution  (Q^).  say  <  1  or  2.  It  is  clear  that  these  [(2)  and  (3)1  are 

conflicting  requirements  and  not  many  investigators  have  taken  full  cognizance 
of  them. "  The  objection  (2)  is  taken  in  failure  to  appreciate  the  logic  of  Knudsen's 
low-pressure  method.  Briefly,  Knudsen's  argument  is  somewhat  as  follows: 

The  molecules  which  at  any  instant  make  up  the  gas  about  a  mean-free-path  A 
from  the  wire  are  those  that  come  from  the  body  of  gas  and  those  that  come  from 
the  wire.  Those  that  come  from  the  body  of  gas  are  brought  up  to  a  cylinder  of 
unit  length  a  radial  distance  ^  from  wire  at  a  rate  about  (2  n  ) );  those  coming 
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from  unit  length  of  wire  do  so  at  the  rate  with  which  they  strike  unit  length  of 
wire,  — -  (2  tt  r^).  Hence,  of  all  the  molecules  which  make  up  the  gas  at  a  dis¬ 
tance  X  from  the  wire  the  fraction  r^/X  come  from  the  wire  carrying  excess  energy 
acquired  at  the  wire.  When  this  fraction  is  small,  say  r^/X^  0.  05,  the  properties 
of  the  gas  at  X  are  practically  uninfluenced  by  the  presence  of  the  wire.  Further¬ 
more,  since  the  molecules  arriving  at  the  wire  have  also  had  their  last  collision 
at  about  X,  the  molecules  arriving  at  the  wire  are  essentially  uninfluenced  by  those 
returning.  Consequently  when  using  cells  of  cylindrical  geometry  the  applicability 
of  free-molecule  heat  conduction  theory  is  tested  by  the  ratio  X/r^  proposed  by 
Knudsen  rather  than  X/r  suggested  in  Reference  (56).  This  is  verified  by  the  expert- 
mental  data  in  Reference  (23).  Consider  the  a  (Kr)  and  a  (A)  values  obtained  by 
Thomas  and  Brown  at  A  T  =  100*^C  in  tube  B,  with  2r  =0.0127  cm  and  2r  =1.95cm. 

1  A 

Q  Q  r\ 

Using  X  (Kr)  =  3.  69  x  10  cm  and  X  (A)  =  4. 71  x  10  cm  at  273  K  and  1  mm  of  Hg 

we  computed  the  Knudsen  numbers  Kn(l)=  X/2r  and  Kn(2)=  X/2r  ,  relative  to 

1  2 

wire  and  tube  diameters,  respectively.  These  are  tabulated  in  Table  II  with  sev¬ 
eral  a  values  given  in  Table  II  of  Reference  (23). 

Table  II 

Accommodation  Coefficients  of  Krypton  and  Argon  on  Platinum  at  T  =  305°K, 

w 

A  T  =  100  according  to  Thomas  and  Brown  with  Corresponding  Knudsen  Numbers. 


(Kr) 

Kn(l) 

Kn(2) 

a  (A) 

Kn(l) 

Kn(2) 

699 

84 

0.55 

.649 

69 

0.45 

686 

74 

0.48 

.648 

46 

0.30 

689 

54 

0.35 

.649 

30 

0.  19 

705 

43 

0.28 

.649 

16 

0.  10 

689 

32 

0.  21 

.640 

12 

0.  064 

693 

26 

0.  17 

703 

21 

0.  14 

665 

11 

0.  073 

634 

6 

0.  038 
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It  is  clear  from  the  table  that  as  long  as  Kn(l)  >  10,  a  calculated  by  the 

low  pressure  theory  is  essentially  independent  of  Kn(l).  The  table  also  shows 

the  low  pressure  theory  is  applicable  (for  the  given  cell  geometry)  with  Kn(2) 

values  about  1/10  the  minimum  suggested  in  Reference  (56).  This  fact  also  renders 

Hartnett's  objection  (3)  less  serious;  since  one  may  work  at  gas  pressures 

about  10  times,  or  more,  greater  than  the  maximxun  he  suggests,  it  is 

possible  to  keep  the  ratio  W  /W  (his  Q^/Q  )  at  a  sensible  value.  For  example, 

with  argon  at  P  =  0.  02639  mm,  and  A  T  =  100°C  (Kn(2)  =  0. 10)  Thomas  and  Brown 

report  W  (filament  area)  =  0. 0147217  watts  and  W  (filament  area)  =  0.  0048545 
S  ^ 

watts;  thus  W  /W  =  3.  0  and  a  =  0.  649.  At  the  lowest  pressure  reported  (0. 00604 

g  V 

mm)  W  would  be  about  0. 23  of  the  above  value  and  assuming  the  same  W  the 
g 

ratio  W  /W  is  about  0, 7.  With  this  ratio  a  was  also  0, 649.  It  is  to  be  noted 
g  V 

that  even  when  working  with  a  gas  such  as  helium  having  an  accommodation 

coefficient  considerably  smaller  than  0.  65  reported  for  argon  in  (^)  one  seldom 

experiences  difficulty  in  finding  a  convenient  pressure  range  within  which  W  is 

g 

a  significant  fraction  of  the  total.power  loss  from  the  filament. 


RESULTS  OF  THEORY  AND  ANALYSIS 

Baule' — ’  reported  the  first  attempt  to  calculate  accommodation  coefficients 
from  known  properties  of  a  gas  and  a  solid.  He  applied  classical  theory  and  used 
a  model  in  which  the  gas  and  surface  atoms  behave  like  elastic  spheres  and  have 
a  Boltzmann  distribution  of  energies.  Baule  supposed  each  gas  atom  to  make 
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a  single  collision  with  a  solid  atom  before  rebounding,  and  applied  the  laws  of 
conservation  of  energy  and  momentum  to  obtain  for  a  the  expression: 


a  = 


2  m  M 


(20) 


(m  +  M) 

where  m  and  M  are  the  masses  of  the  gas  and  surface  atoms.  In  case  the  surface 

is  rough  and  multiple  collisions  occur,  say  s  in  number,  then  according  to  Baule 

the  apparent  (or  observed)  accommodation  coefficient,  a  is  given  fcy: 

-nS  s 

2  2 


a^  =  l- 


m  +  M 


(21) 


,125) 


(m  +  M) 

Equation  (20)  is  derived  in  Reference  (T_),  Equations  (21)  and  (20)  combined 
result  in  an  expression  identical  in  form  with  the  formula  given  by  Roberts  to 
account  for  the  effect  of  surface  roughness  on  measured  a ;  viz: 

1/s 


0,  =  1  -  ( 1  -  a  ) 
s 


(22) 


Brown^  — ^  compared  the  variation  of  the  ratio  a  (Ne)/  a  (He)  with  M  predicted 

by  Equation  (21),  with  experimental  data  in  somewhat  the  following  way:  Equation 

s  ^ S"*!  y  2 

(21)  may  be  written  as  a  =  s  a  - — ^  a  +  .  .  .  ;  and  since  a(He)  and  a  (Ne) 

s  *  • 

on  clean  metal  surfaces  are  considerably  less  than  unity  the  Equation  (21)  can 
be  reduced  to: 


=  s  =  2sm  M/  (m  +  M)^  (23) 

s 

Figure  (7)  according  to  Reference  (10)  shows  a  plot  of  ^  (Ne)  /  a  (He)  vs.  M 

computed  from  experimental  data  and  from  Equation  (23)  in  which 

taken  as  unity  on  the  assumption  that  s  depends  only  on  the  structure  of  the  surface. 
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( 58 ) 

Hurlbut' — '  outlined  objections  to  the  Baule  model  which  are,  in  part,  as 
follows:  The  choice  of  a  classical  model  is  questionable,  and  the  assumption 
that  each  solid  atom  is  in  independent  elastic  suspension  is  not  realistic.  Exper¬ 
imental  results  suggest  that  a  surface  possesses  greater  effective  mass  in  colli¬ 
sion  than  possessed  by  any  one  surface  atom.  The  assumption  that  the  time  of 
collision  is  infinitely  short  is  also  doubtful,  especially  at  high  interaction  energies 
where  an  incident  particle  may  penetrate  below  the  surface  of  the  lattice.  In 
addition  to  these  is  the  empirical  fact  that  experimental  values  of  a  (He  and  Ne 
on  W,  for  example)  lead  to  fractional  numbers  of  collisions  with  die  surface, 
(Equation  (21))  a  result  in  conflict  with  reality. 

Modem  theory  begins  with  Zener^~^  and  Jackson^  — ^  et.  al.  *  (1^)^ 

The  model  is  of  gas  atoms  which  have  a  one  dimensional  distribution  of  energies. 
The  lines  of  motion  of  solid  atoms  and  incident  gas  particles  are  assumed  normal 
to  the  surface.  The  theoretical  approach,  in  general,  is  to  consider  the  solid  atoms 
as  oscillators  which  behave  in  the  same  way  as  atoms  in  the  interior  of  the  solid 
and  compute  the  probability  that  a  solid  atom  in  vibrational  state  i  struck  by  a 
gas  atom  with  energy  E  will  undergo  a  transition  to  state  i' ,  the  gas  atom  depart¬ 
ing  with  energy  E  +  (i-i*)hv.  Jackson  and  Mott^~^  considered  the  energy 
exchange  to  be  nearly  always  limited  to  a  single  quantum  transition  of  the  oscil¬ 
lator.  The  assumption  is  that  when  the  interaction  energy  of  a  gas  atom  with 
the  solid  is  expanded  in  powers  of  displacements  of  lattice  atoms  from  their 
equilibrium  positions  only  linear  terms  are  important.  Thus  in  the  final  formulas 
for  a,  ji  -  i*  1  is  always  taken  as  unity.  Strachan^— ^  has  justified  this  assumption 
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by  showing  that  the  probability  of  a  solid  changing  a  given  amount  of  energy  by  a 
one-quantum  jump  is  much  larger  than  the  probability  of  a  change  by  a  two  quan¬ 
tum  jump. 

One  important  difference  among  the  several  approaches  to  theory  is  in  the 

treatment  of  the  solid.  In  the  earlier  theory  the  distribution  of  energy  amongst  the 

atoms  of  the  solid  is  taken  to  be  that  of  an  assembly  of  oscillators  all  of  the  same 

frequency  v  as  in  Einstein's  specific  heat  theory.  In  later  theory  account  is  taken 

of  the  detailed  distribution  of  normal  modes  within  the  crystal,  i.e. ,  the  solid  is 

treated  in  the  same  way  as  in  the  Debye  specific  heat  theory.  The  accommodation 

coefficient,  as  defined  by  Equation  (3),  is  calculated  from  the  transition  probability 

by  averaging  over  the  energies  of  the  gas  atoms  and  over  the  vibrational  fre- 

(12) 

quency  distribution  of  the  sblid  atoms.  Jackson  and  Howarth' — 'indicate  how  this 
average  may  be  performed. 

Another  important  difference  among  theoretical  treatments  is  in  the  assump¬ 
tions  regarding  the  interaction  between  gas  atoms  and  solid  surface.  Jackson 
and  Mott^— ^  considered  two  types  of  interaction;  In  one  a  one-dimensional 

repulsive  exponential  field  was  assumed,  with  potential  energy  at  an  interatomic 

-ay 

distance  y  of  the  form  V  =  C  e  ,  where  C  and  a  are  constants;  in  the  other  a 

rigid  elastic  sphere  model  was  treated;  taking  y  as  the  distance  of  closest 

o 

approach  of  gas  and  solid  atoms  V  was  set  at  zero  for  y  >  y  and  at  infinity 

o 

for  y<y^  •  In  the  latter  case  the  expression  for  a  does  not  contain  adjustable 
constants;  it  results  in  a  ratio  ^  (Ne)/  (He)  about  5  on  the  same  metal  surface 
while  Roberts'  values  give  1.  2  and  more  recent  values  about  2.3.  Calculations 
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of  a  values  as  a  function  of  temperature  according  to  the  theory  (with  a  =  9  x 
8  “1 

10  cm  )  result  in  satisfactory  agreement  with  the  temperature  dependence 
exhibited  by  Roberts'  a  (He  -  W)  values  These  three  clean  surface  a 

values  were  the  only  ones  available  at  the  time  the  Jackson-Mott  theory  was 
developed. 

-ay 

The  potential  function  Ce  implies  repulsion  at  all  distances.  It  is 
known,  however,  that  there  is  an  attractive  field  between  a  solid  and  a  gas  (the 
fact  that  adsorption  can  take  place  is  evidence  that  such  a  field  exists).  The 
attractive  field  is  considerably  greater  for  neon  than  it  is  for  helium,  and  the 
theory  described  thus  far  is  not  applicable  to  experimental  results  with  neon. 

r  61  ^ 

Devonshire'  — '  in  one  of  a  series  of  papers  by  Lennard-Jones  and  his  co-workers 
attempted  to  account  for  the  attractive  as  well  as  the  repulsive  field  by  assuming 
an  interaction  potential  in  the  form  of  a  Morse  potential  function.  In  this  form 
the  interaction  potential  is  attractive  at  large  distances  and  repulsive  at  small 
ones  with  a  minimum  (potential  well)  in  between.  Devonshire  did  not  treat  the 
possibility  that  a  gas  atom  may  be  entrapped  in  vibrational  states  in  the  potential 
well.  However,  the  formulas  for  a  which  Devonshire  derived  are  more  general 
than  previous  ones  ^  — ^  ^  and  contain  them  as  special  cases.  Comparison 

of  theoretical  computations  with  the  temperature  dependence  exhibited  by  Robert's 
a  (He  -  W)  and  a  (Ne  -  W)  values  gives  satisfactory  agreement.  Brown^~^  showed 
how  Devonshire's  formula  may  be  used  to  obtain  agreement  with  more  recent  experi- 
mental  values  of  a  (He)  and  a  (Ne)  on  tungsten'^ —  ,  iron' —  ,  aluminum'  —  ^  and 
platinum^  — ^  However,  in  no  case  can  a  be  calculated  a  priori  from  theory  since 
agreement  with  experimental  data  is  obtained  only  after  judiciously  adjusting 
several  constants  in  the  final  formulas. 
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Devonshire  also  obtained  the  result  that  a  is  proportional  to  1/ 6  where  0 
is  the  characteristic  temperature  of  the  solid.  Brown^— ^  showed  that  upon 
making  assumptions  regarding  the  potential  field  which  are  necessary  to  achieve 
agreement  between  the  Devonshire  theory  and  experiment  it  is  found  that  a  is 

3 

more  nearly  proportional  to  1  /M  0  where  M  is  the  molecular  weight  of  the  solid. 
Brown  also  observed  that  Devonshire  did  not  consider  the  difference  that  can 
exist  between  the  attractive  and  repulsive  parts  of  the  Morse  potential  function. 
Devonshire  considered  both  parts  to  vary  with  the  motion  of  the  solid  atom  where¬ 
as,  as  pointed  out  by  Jackson^— \  the  attractive  part,  which  particularly  at 
large  distances  depends  on  the  average  properties  of  a  considerable  volume  of 
the  solid  in  the  neighborhood  of  the  atom  struck,  should  be  stationary.  Accordingly, 
Brown  showed  how  Devonshire's  treatment  can  be  modified  so  as  to  disengage 
the  attractive  term  in  the  Morse  function  from  the  oscillating  solid  atom. 

Devonshire  examined  the  effect  on  a  of  using  a  one-dimensional  model 
rnd  surmised  that  the  result  is  an  a  value  which  is  too  large  by  a  small  amount. 

The  theory  which  has  been  described  applies  only  to  metallic  surfaces 
which  are  reasonably  free  of  adsorbed  gases.  Hence,  it  is  only  applicable  to 
interactions  between  a  metal  surface  and  a  gas  which  does  not  adsorb  on  that 
surface.  Michels'  — '  attempted  to  relate  a  on  a  clean  surface  to  a  on  a  con¬ 
taminated  surface.  However,  to  date,  there  is  no  known  quantitative  functional 
relation  between  n  and  the  frac*^ion  of  the  surface  which  is  covered,  a  relation 
which  may  be  nonlinear. 

We  have  thus  far,  in  the  discussion,  failed  to  note  that  the  accommodation 
coefficient  for  internal  energy  n  .  may  be  different  than  that  for  translational 
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/  9 ) 

energy  a  ,  and  assumed  with  Knudsen  —  that  a ,  and  a.  have  a  common  value  a . 
t  r  1 

According  to  Knudsen  the  effect  of  internal  accommodation  can  be  taken  into 
account  by  the  relation: 


W 


a  +  T  f  a .  = 
t  4 


(24) 


^  [  P/y2nMRTg]  2R(^T^-T^j 

where  the  r.h.  s.  is  obtained  from  Equations  (9)  and  (12)  and  the  factor  f  depends 
on  the  molar  heat  capacities  and  at  constant  pressure  and  volume,  respec¬ 


tively,  and  is  given  by: 
f  = 


5/3  -  C  /C 

£_  V 


C  /C  -  1 
P  V 


(25) 


When  a  ^  =  a  ^ ,  then  the  expression 


a  = 


yR/2nMTg  (C^R  +  1/2)  (T^  -  T^) 


(27) 


is  obtained  which  is  identical  with  the  formula  for  a  that  results  from  combining 
Equations  (9)  and  (13). 

Knudsen's  assumption  that  a  ^  =  a .  =  a  is  supported  by  results  of  his  own 
experimental  observations;  however,  there  is  also  some  contradictory  evidence: 
Rowley  and  Bonhoeffer^~^  measured  accommodation  coefficients  of  hydrogen 
on  platinum  and  compared  the  accommodation  coefficient  of  normal  hydrogen 
with  that  of  para-hydrogen  at  temperatures  where  p  -  H  has  only  translational 

Ci 

energy  but  n  -  H  still  has  rotational  energy  in  its  ortho  component.  Measure- 

ments  resulted  in  a  ^  =  0.  43  and  a  ^  =  0.  26  at  140°K  and  a  ^  =  0.  44,  a  ^  =  0. 18 
t  rot  t  rot 

at  170°K,  Others  (65)  -  {§§_)  and  (^),  (70),  (n.  v. )  have  questioned  the  assump¬ 
tion  that  ^  j  proposed  various  analytical  and  experimental  methods  for 
determining  a  for  polyatomic  gases. 
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( 71 ) 

Jackson  and  Howarth  —  attempted  to  modify  the  collision  theory,  pre¬ 
viously  described,  to  include  the  effect  of  rotation  of  diatomic  molecules.  They 
considered  only  repulsive  interactions  and  two  limiting  cases:  In  one  the  gas 
molecule  is  treated  as  an  oscillator,  in  the  other  as  a  plane  rotator.  In  the 
first  case  the  energy  exchange  between  the  rotations  of  the  diatomic  molecule 
and  the  vibrations  of  the  solid  is  negligible;  in  the  second,  assuming  the  same 
repulsive  field  for  the  solid,  the  effect  of  rotation  is  small  for  Hydrogen  and 
larger  for  Oxygen. 

Other  theoretical  work  includes  (22)  -  { 75 )  (n.  v. )  and  a  review  paper  by 
Herzfeld^  — 

The  energy  exchange  between  a  gas  and  a  solid  which  takes  place  when  a 

gas  particle  becomes  adsorbed  on  the  surface  first  and  re-evaporated  at  a  later 

time  has  also  been  treated  theoretically — '  ' — The  distinction  between 

this  type  of  an  interaction  and  the  fype  considered  previously  viz:  an  inelastic 

collision  followed  "immediately"  by  reflection,  is  not  always  well  defined.  The 

(77) 

assumption  sometimes  made  —  that  all  impinging  gas  atoms  are  temporarily 
adsorbed  has  not  been  generally  demonsti'ated,  and  may  be  questioned  in  the 
case  of  helium,  for  example,  interacting  with  a  clean  metal  surface.  The  fact 
that  the  accommodation  coefficient  of  helium  on  a  clean  surface  appears  to  be 
unaltered  by  considerable  variation  in  A  T  at  constant  Tg  (see  above)  suggests 
that  at  least  for  helium  adsorption  does  not  play  a  major  role  in  the  energy  ex¬ 
change  process. 
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CONCLUDING  REMARKS 


Accommodation  coefficient  values  which  can  be  confidently  applied  to  a 
solution  of  satellite-drag  and  heat -transfer  (80)  problems  must  be  obtained 
from  experiments  with  gas  molecules  whose  velocities  match  those  of  satellites. 

To  our  knowledge,  the  problem  of  generating  in  the  laboratory  a  nearly  mono- 
energetic  beam  of  molecules  in  the  range  1  -  10  eV  with  adequate  flux  has  not 
yet  been  solved.  Experimental  facilities  such  as  shock  and  arc  tunnels  which 
may  produce  molecular  flows  of  nearly  adequate  velocities  are  not  readily 
adaptable  to  a  measurements.  For  these  measurements  the  transient  flow 
obtained  in  a  shock  tunnel  is  extremely  inconvenient,  while  the  composition  and 
energy  distribution  of  the  gas  obtained  in  an  arc  tunnel  is  not  known  with  suffi¬ 
cient  accuracy.  With  both  facilities  the  problem  of  characterizing  the  composi¬ 
tion  of  the  surface  is  formidable.  One  may  question  the  need  for  control  of 
surface  composition  in  experiments  designed  for  satellite  applications.  How¬ 
ever,  the  composition  of  a  satellite  surface  in  orbit,  and  hence  the  accommoda¬ 
tion  coefficient,  will  most  probably  differ  appreciably  from  that  of  an  untreated 

-7 

surface  of  the  same  material  exposed  to  a  laboratory  vacuum  of,  say  only,  10  mm 
of  Hg.  Furthermore,  in  orbit,  satellite  surface  composition  may  change  with 
changing  surface  temperature  and  ambient  gas  pressnrt-,  resulting  in  corres¬ 
ponding  changes  in  a  .  Thus  it  may  well  be  necessary  to  i  .rntrol  in  the  laboratory, 
the  composition  of  the  surface  in  order  to  obtain  useful  data  for  satellite  appli¬ 
cations.  Such  control  will  certainly  be  necessary  if  the  data  are  to  be  used  to 
test  theory. 
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